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Abstract. The presence of surfactants has a pronounced effect on the surface 
tension and, hence, on the stress balance at the phase separating interface of 
two-phase flows. The transport of momentum induced by the local variations 
of the capillary forces are known as Marangoni effects. Here we study a model, 
which assumes the surfactant to be soluble in one of the adjacent bulk phases 
and which represents a generalization of the two-phase Navier-Stokes equa- 
tions. Based on maximal Lp-regularity results for suitable linearizations we 
obtain local well-posedness of this model. We employ recent results from the 
Lp-theory of two-phase flows without surfactant. 



Introduction 

We consider a free boundary problem, which describes the isothermal flow of two 
immiscible, incompressible Newtonian fluids. To be precise, we assume the fluids to 
have constant densities p± > and constant viscosities ri± > Q. They occupy the 
bounded domain Q. C M", where the dispersed phase is located in C Vl and 

separated from the continuous phase, which is located in ^+{t) C f2, by a sharp 
interface T{t). The phase configuration is completely determined by the position 
of the interface, which, however, is time dependent and has to be resolved as part 
of the problem. The evolution of the velocity field u and the pressure p is governed 
by the Navier-Stokes equations 

pD-^u-riAu + Vp ^ pf, t>0,xen±{t), 

divu = 0, t > 0, X e fl±{t), 

Man = 0, t>0, xedn, 

u{0) = uo, X e ri±(o), 

where we denote by D" = dt + u ■ ^ the material derivative w. r.t. u and [•]dn 
denotes the trace of a quantity defined in ^+{t) on the boundary 9r2. As above 
it will be convenient to drop the phase subscripts and to write p, rj, ... instead of 
P± I 'yi I ■ • • J whenever there is no danger of confusion. 

We assume, that no phase transitions and no interfacial slip occur, which implies 
the velocity field to be continuous across the interface. In this case the normal 
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velocity Vr of the interface equals the normal component of the velocity field. In 
summary, the initial phase configuration and the transmission conditions for mass 
and momentum 

Mr(t) =0, t>o, xe r{t), 

-MVu + VvJ) - i3]r(t) i^nt) = divr(t) { (TPr(t) }, t > 0, x € r{t), 

(lb) 

= Mr(t) • T^m, t>0, xG T{t), 

r(o) = To 

completely determine the evolution of the interface. Here, i^r(t) denotes the normal 
field on r(<) pointing from Q-{t) into n+{t). Moreover, Pr{t) = 1 — ^'r(t) ® 2^r(t) 
denotes the projection onto the tangent space of T{t) and [ • ]r(() denotes the trace 
on T{t) of a quantity defined in Q±{t). Finally, |-lr(t) denotes the jump of a 
quantity defined in fl±{t) across T{t), i.e. 

Mnt) [t, x) = 

J^^{ <l>{t, X + efr(t)(i, a;)) - 4>{t, x - efr(t)(i, a;)) }, t> 0, x e T{t). 

If a surface active agent - surfactant for short - is present, then the surface 
tension fj > depends on the surface specific concentration of this surfactant via a 
so-called equation of state 

(Ic) a = a(cr). 

Hence, we have 

divr(t) { a{cr)Pr{t) } = f^(cr) Kr(t) 2^(t) + cr'(cr)Vr(t) cr, 

where K-p(i) denotes the sum of the principal curvatures of the interface. Note that 
the surfacant is adsorbed at the interface r(<) in this case and we denote by cr its 
surface specific concentration. 

Moreover, we assume such a surfactant to be soluble in the continuous phase 
^+(t) and denote by c its volume specific concentration. Based on Ficks law of 
diffusion the evolution of the surfactant is governed by 

D^c - dAc = 0, t > 0, ;r e 

a([c]r(i)) = Cr, t > 0, x e r(i), 

(Id) D^cr + Cr divr(t) u - dr Ar(t)Cr = dd+c, t> 0, x e T{t), 

= t>o, xe on, 

c(o) = Co, xen+{0), 

where the normal derivatives have to be understood as = t[V • ^]r(t) ^'r(t) resp. 

= [\/ ■ ^]an i^dQ- Note that cr is defined only on the graph gr{T{t)). However, due 
to the kinematic condition Vr = [u]r{t) " ^(t) its material derivative is well-defined 
as 

D^crit, x) := :^cr(s, x"(s \t, x)) , t>0, xe T{t), 
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where x"( ' I denotes the characteristic curve of a particle, which at time t > 
is located at a; e 57 and which is advected by the velocity field u, i. e. 

r{s\t,x)^u{s,x"{s\t,x)), |s-t|<e, x''{t\t,x)^x. 

In fact, X £ r{t) implies x"('S 1^1^) ^ r(s), if |s — <| < e. This is the reason why 
we prefer the notation employing the material derivative even for the momentum 
balance in (llap and the surfactant diffusion equation in the bulk phase in (jldp . 

1. An Analytic Approach 

The two-phase Navier-Stokes equations with surface tension, i. e. (ITalllb|) with 
constant cr > 0, have already been intensively studied. For the most recent results 
concerning local-in-time well-posedness as well as the qualitative behaviour of so- 
lutions we refer to jU [5] and the references therein. The two-phase Navier-Stokes 
equations with soluble surfactant ([T|) have been studied in a prototype geometry, 
where the interface is almost flat, see 0. The first aim of the present paper is to 
transfer this result to a general geometry, i. e. we will prove 

Theorem 1.1. Let fl C M" be a bounded domain with boundary of class and 
let p > n + 2. Let p±, r]±, d, dr > and let a, a e C^~(R+, M+) with a' > 0. 
Suppose 

no e w^'^/^in \ To), To e w^-'/p, CO G w^-'^^in+io), M+) 

are subject to the regularity and compatibility conditions 

divuo^O inri\ro, ImqIfo = 0, [wo]ao = 0, 

--PFo[»y(Vw + VuT)]r„ = Vro (o-o ao [co]fo), 

[co]fo e W^-^/^iTo), d,co = on 8^ 

and f e Lp(R+, Lp{D,, R")). Then there exists oq — ao{uo, To, cq) > 0, such that 
the two-phase Navier-Stokes equations with soluble surfactant ([IJ) admit a unique 
local strong solution (u, p, T, c, cr) with c, cr > in (0, ao). 

If in addition a, a € C"(R+, R+) and f — 0, then the above solution is a 
classical solution. The graph 

grr= U {t}xT{t) 

0<t<ao 

is a real analytic manifold and with 

I5 = {{t,x)e (0, ao)xn : x(^ T{t) }, 15+ ^ { {t, x) € (0, ao) x fl : x £ n+{t) } 
the functions {u, p) : 13 — > R"+-'^ and c : 0+ — > R+ are real analytic. 

The proof of Theorem 11.11 will be carried out in Sections [3] and SI However, 
besides the local well-posedness of model (H]) it is worthwhile to study its dynam- 
ics, i. e. to characterize the equilibria and to study the qualitative behaviour of the 
solutions. The state of the system is described in terms of the velocity, the position 
of the interface and the surfactant concentration in the bulk phase. To determine 
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the corresponding phase manifold, we denote by A4H^{il) the manifold of all com- 
pact C^-hypersurfaces in fl, which becomes a metric space when equipped with the 
Hausdorff metric 



dist / 



)(E, r) :=dist(iV2S, N'^T) 



s, r e Mn^{n). 



Here, N'^Y, denotes the second normal bundle of a hypersurface E G AiH^{il). As 
has been noted in [H Section 5], every hypersurface E € MH^{^) may be described 
by a level-set function ip^: G C^(n), i. e. S = '^^^{0), and, of course, S is of class 
Wp, if and only if ip^ G Wp{V,). With this notation, the phase manifold of the 
dynamical system ([!]) is given by 

{v, E, x) e C(n)" X X C'(f^) 



(2) Spin) 



divv = in \ E, Hs = 0, [v]9n = 0, 



domx = n\ int(E), x e Wp "'''{n \ int(E)), 



,2-2/p 



on 917 
(7'(a([x]s))VEa([x]E) J 



where we denote by 

c{n) 



{ (f) G C([/) : [/ C open, 



C/C f7} 



the set of functions, which are defined and continuous on certain subsets of fJ, 
and dom (p denotes the domain of definition of such a function </>. Moreover, for 
E G we denote by int(E) the interior of E, i. e. the union of all connected 

components of £7 \ E, which are not in contact with dVl. Since the solutions to ([T|), 
which exist thanks to Theorem 11.11 preserve the compatibility conditions of the 
phase manifold, we may prove 

Theorem 1.2. Let Q, C R" he a hounded domain with boundary of class and 
let p > n + 2. Let p±, ri±, d, dr > and let a, a e C^~{R+, R+) with a' > 0. 
Then the strong solutions to {IJ) generate a local semiflow in the phase manifold 
Sp{n). Each of these solutions exists on a maximal time interval [0, a*). 

The proof of Theorem 11.21 will be carried out in Section [5] Note that neither 
the pressure nor the surfactant concentration on the interface explicitly appear in 
the definition of the phase manifold Sp{VL). In fact, the surfactant concentration 
on the interface is well-defined for every state (u, E, x) G Sp{n) as xr = a([x]s)- 
Moreover, the pressure may be reconstructed for every semiflow 

[u, r, c) G BUC{[Q, ao), Sp{n)), 

which is induced by a strong solution to ([Ij, i.e. which in particular satisfies 
9tu, Am G ip((0, ao) x fi)". Indeed, in this case we have 

(u|V0)n = -(div?/|0)o-(Mr-i'r|0)r=O, G iJ^M^!), 



ON TWO-PHASE FLOWS WITH SOLUBLE SURFACTANT 



5 



which implies {dtu\W(p)n = for ah (p £ Hp,{il), where l/p+ 1/p' = 1. Hence, 
the momentum balance in (|lap and the momentum transmission condition in (|lb[) 
yield 

(V7r|V0)n = {^Au~u-Wu\Wcb)n, cf, e Hl,{n), 

Mr(t) = cr(cr) Kr(t) + Ivi^u + Vu^)jr(t) '^r(t) ■ ^'r(t), * > 0, x e r{t) 

for the modified pressure tt — p/ p, where, of course, cr = a([c]r). Thanks to [6l 
Theorem 8.5], this problem has a unique solution tt G ip((0, oq) x such that 
7r(i) e \ r(t)) for < t < ao- 

As has been worked out in [1], the equilibria are closely related to the energy 
functional 

(3a) mu,T, c)^\\\^pu{t)\\l^^^^^^)+ J 0(c)dH"+ J 0r(a([c]r(t)))dH"-\ 

n+(t) r(t) 

which is composed of the kinetic energy and the free energies in the bulk phase and 
on the interface. The latter are obtained via the energy densities (f> and (pr, which 
may be obtained as 

s s 

(3b) 4>is) = j pT{a{r))Ar, (/)r(s) = cr(s) + s/ir(s), -A*r(s) = Ja'{r)/rdr 



based on the chemical potential pr- With these definitions T, Theorem 3.1] states 
the following. 

Theorem 1.3. Let fl C M" be a bounded domain with boundary of class C^^ and 
letp > n + 2. Let p±, ri±, d, dr > and let cr, a G C'^^(M+, M+) with —a', a' > 0. 
Let the energy functional $ and the energy densities (f), (f>r be defined by Then 
we have 

L The energy equality 

9*$ + 277 j |i:'pd-H" + d j cj)"{c)\\7c\^dK' + dr J p'^icr) \\/rcr\^ dV.''-^ = 
o\r(i) n+(t) r(t) 

is valid for smooth solutions to (Qp. 

2. The equilibria of {Ip are zero velocities, constant pressures in the connected 
components of the two phases, constant surfactant concentrations in fl^ and 
on r and the dispersed phase is a union of non-intersecting balls. 

3. The energy functional ^ is a strict Lyapunov functional. 

4. The critical points of the energy functional $ for prescribed volumes of the con- 
nected components of the dispersed phase and prescribed total surfactant mass 
are precisely the equilibria of the system. □ 

Note that (j3bp implies (f>"{s) = —a'{a{s))a'{s)/a{s) as well as (/'"(s) = —cf'{s)/s 
for s > Q, i.e. (j) and (pr are convex, provided the surface tension a is monotonically 
decreasing and the adsorption isotherm a is monotonically increasing, which are 
natural assumptions. 
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2. Linearization 

2.1. The Hanzawa Transformation. The proof of Theorem ll.ll is based on max- 
imal regularity results for suitable linearizations of ([T]) and fixed point arguments. 
To obtain a linearization in a fixed domain note that any hypersurface S G M'H^{fl) 
forms the boundary of the union of finitely many domains contained in f2, i. e. 
int(S]) C is an open set with compact boundary of class and outer unit nor- 
mal field v^. Moreover, every such hypersurface admits a tube contained in 17, i. e. 
there exists e > 0, such that the tube mapping 

r| : (-e, e) X S — > Q, T{r, x) = x + ri^j:{x), |r| < e, x e E 

is a C^-diffeomorphism from (— e, e) x S onto R{T^) C fi. The supremum of all 
such e > is called the tube size of S and denoted by t{Yj). Note that the tube 
size is bounded from above by the radii of the interior and exterior ball condition 
of int(S) as well as by the reciprocals of the moduli of the principal curvatures 
Ki, . . . , Kn-i, which are the eigenvalues of the curvature tensor Ls = — Vs i^s. In 
particular, we have 

r(E) < min | l/|Ki(x)| : a: G S, « = 1, . . . , n — 1 }. 

Now, given any initial interface Fq G Wp and any (5 > 0, it is well-known 
that there exists a reference hypersurface S G A4H^{fl), which may be choosen to 
be real analytic, such that 

dist^«2(n)(E, To) <(5 and To e R{T^^^) 

for some < e < t(E), cf. [H Section 2]. Hence, Fq may be parametrized over E 
via a function 70 G W^~^^^(E) as 

To = {x + jo{x)iy^{x) : a; G E} 

and we expect the interface T{t) to stay in the tube R(T^^^) and to be parametrized 
by a function j(t, ■ ) G Wp ^^''(E) as well, at least for small times < t < a. Thus, 
by choosing a cut-off function x G C'o°(I'^; [0, 1]) with 

X(r) = l, |r|<l/3, xW=0, |r| > 2/3, Ix'U < 4 

we obtain a family of C^-diffcomorphisms, the Hanzawa transformation, 

e : [0, a) X — >n, 
e(t, x) = x + x{dT.{x)ie) "fit, n^{x)) z/s(nE(2:)), < t < a, x € n, 

which are of class Wp if 7 G BUC{[0, a), Wp ^/^(E)) and E is choosen to be 
sufficiently smooth, and which transform the system ([T]) to a fixed domain, since 

n^{t) = int(r(i)) = e(t, int(E)), n+{t) = SI \ iHt(r(t)) = e{t, f7\5t(E)) 

and, of course, r(<) = 9(t, E) for all < i < a. Note, that we denote by 
T^" = ds X Hs : R{Tl) — > (-£, e) x E 
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the inverse of T^, i.e. : R{T-^) — > {~£, e) is nothing but the signed distance 
from S and : R{T^) — !• E is the metric projection onto S. Also note, that 
the transformation Q{t, ■ ) is weh-defined for aU x S fi, even if o?s and Us are only 
defined in R{T^), since the cut-off function x{ • /e) vanishes outside (— 2e/3, 2e/3) 
and 9(t, • ) is to be understood as the identity outside R{t'^^^). 

For the following computations it will be convenient, to have an alternative 
representation of the Hanzawa transformation at hand. With < t < a and x e S 
fixed we define 

0{-\t, x) : {-e, e) — > (-£, e), e{r ; x) -.^ r + x{r /e)-i{t, x), \r\ < e 
and note that 

(4a) e(t, x)^x + { e{ds{x) ; t, Us{x)) ~ ds{x) } lysi^six)), < t < a, x e fl. 

Thus, for fixed < i < a and a; e S, the mapping 9{- ; t, x) : (— e, e) — !• (— e, e) is 
a bijection and we denote its inverse by ^( • ; t, x) : (— e, e) — > (— e, e). Note that 
9{r ; t, x) = 9{r ■ t, x) = r, if \r\ > 2e/3 and, therefore, we may extend 0{ - ; t, x) 
and 6{- ; t, x) as the identity to R \ (— e, e) and obtain a pair of inverse bijections. 
Since S is compact, the signed distance extends to a functional d-^ '■ fl — > M via 

=Fdists(x), a; e int(I]) U S resp. a; e ext(E) U S, 

where we denote by ext(I]) ;= 51 \ int(I]) the continuous phase in the transformed 
setting. Thus, the spatial inverse Q(t, • ) : 57 — !• 57 of Q{t, • ) : 57 — > 57 is given as 

(4b) e(t, a;) = a; -f { e{dsix) ; t, Us{x)) - ds{x) } i^s(ns(a;)), < t < a, a; e 57. 

Note that in both representations (|4]) the projection onto S may be interpreted as 
the multi- valued metric projection Us : 57 — > 2^ \ { }, since 

< t < 0, y G S, 

0(ds(x) ;t,y)~ ds(x) = 0{dsix) ; t, y) - ds{x) = 0, 

a; e 57, \dsix)\ > 2e/3 

and lis : R{T§^) — ^ ^ is well-defined as a single-valued metric projection. 

Finally note that the function 0{ - ; t, x) : K — >■ M is monotonically increasing 
for all < t < a and a; G S, provided that \\^{t, ■ )\\bc{S) < ^/4. This shows that 
Q(t, • ) : 57 — > 57 is indeed invertible and we have 

0{r \t, x) ^r- x(r/£ ; t, x) -/{t, x), r G R, < t < a, a; G E 

for some function x( ■ ; 2;) : M — s- R with xi^ t, = 0, if \r\ > 2/3. Moreover, 
we have 

Ve(i, a;) = 1 + lx'{d^{x)/eh{t, ns(a:))(l - Ps(ns(a;))) 

+ x{dsix)/eMdsix)){S/sl){t, Us{x)) ® M^si^)) 

- x{d^{x)/sh{t, ns(x))N(ds(x))Ls(ns(x)) 

for < i < a, a; G 57, which shows that V8(i, • ) is boundedly invertible and &{t, • ) 
is indeed a C^-diffeomorphism, provided \\j{t, • )||_bc(i;) ^^id ||Vs7(i, • )||_bc(i;) 
sufficiently small. Here N( • ) = (1 — • i^)^^ is defined by (|6]). 
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2.2. The Transformed Equations. For the transformed quantities 

u{t, ■ ) := u{t, ■ ) o e{t, ■ ), p{t, ■ ) := p{t, ■ ) o e{t, ■ ) 
c{t, ■ ) := c{t, ■ ) o [e{t, ■ )]ext(E) 

we obtain transformed differential operators via 



£(7) 4>{t, x) 



= {D'^(j)}{t,e{t,x)), 0<i<a, a;eO\S, 
= {V(^}(i, G(t, a;)), 0<t<o, a;eO\I], 
= { A(^} (f, e(i, a;)), Q <t < a, x &Q.\Y,, 



where ^(i, • ) = (t>{t, ■ ) o Q{t, • ) : O \ S 
scalar valued function (/>, as well as 



denotes the transform of a generic 



D(7)$(i, x) := {div$} (i, Q{t, x)), < t < a, a;eO\S, 



where • ) = • ) o Q{t, • ) : \ S 
vector field Moreover, we set 



I" denotes the transform of a generic 



M^{j) Mt, x) 
£r(7) 4'^{t, x) 



= {D^(l)r}it, e{t, x)), 0<i<a, xeS, 
= {Vr(t)0r}(i, e(t, x)), < i < a, x € S, 
= {Ar(t)?ir}(t, e(t, x)), 0<t<a, xeS, 



where (/isCi, • ) = (f>r{t, ■) o Q(t, • ) : E — > M denotes the transform of a generic 
scalar valued function <pr, and, analogously. 



I?r(7) ^E(i, a;) := { divr(t) $r } {t, &{t, x)), < t < a, x G E, 



where ^^{t, ■) = ^r{t, ■) o Q[t, ■) : E — > R" denotes the transform of a generic 
vector field $r- Finally, we transform the geometric quantities 



Ml) (*> x) 

Kr(7) {t, x) 
Vr(7) (i, a;) 



:= fr(t)(i, a;)), 
:= Kr(t)(i, ©(i, a;)), 
:= yr(i, e(i, a;)), 



< t < a, a; e S, 
< t < a, a; e S, 
0<t<a, a;eS 
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and arrive at 

p7W"(7)?2 - vC{"f)u + g{"f)p = pf in (0, a) X \ E, 

V{-f)u = in (0, a) X f7\ S, 

A^"(7)c-d/:(7)c = in (0, a) x ext(S), 

|m]e = 0, Vr{j) = u ■ i/r(7) on (0, a) x E, 
-[77(^(7)^ + ^(7)^"^) -p]Ei^r(7) 

(5) 

-I?r(7){ 'T(cs)Pr(7) } = on (0, a) x E, 

"([c]e) = cs on (0, a) x E, 

A^"(7)cE + csl?r(7)" - rfr/:r(7)cs = d^^c on (0, a) x E, 

= 0, 9i,c = on (0, a) x 9fi, 

u(0) = uo in ri, 7(0) = 7o on E, c(0) = cq in ext(E), 

where d1 on denotes the normal derivative on E = c)ext(E) n E, which should be 
interpreted as part of the boundary of ext(E) in this case, i. e. the normal points 
from ext(E) into int(E). Furthermore, the operator Pr(7) = 1 ^ ^r(7) ® i^r(7) 
denotes the transformed projection ^r(.) and 

CE(t, •) :=cr(i, ■)oQ(t, •) 

denotes the transformed surface specific surfactant concentration, being a function 
ce : (0, a) X E — M+, which explains the change in the subscript. 

Note that problem ([5]) may be interpreted as an evolution equation for the trans- 
formed unkowns u and c in the fixed domain 17 \ E resp. ext(E). The pressure p 
plays the role of a Lagrangian parameter introduced due to the divergence con- 
straint. Moreover, the velocity u is subject to a dynamic transmission condition 
across E, which introduces the parametrization 7 and the corresponding evolution 
equation on E. Finally, the surfactant concentration c is subject to a dynamic 
boundary condition on E, which introduces the surface concentration cs of ad- 
sorbed surfactant and the corresponding evolution equation on E. Indeed, we have 

vv{-i) ^ Ai(7)(f^E - N(7)Ve7), 

6 

/.(7) = (1 + |N(7)Vs7p)-'/^ N(7) - (1 - lL^)-\ 

where L^, — Vsi^s denotes the curvature tensor of E, and the kinematic condition 
reads 

da = At(7)"^9t7(z/s • t'r(7)) = Kl)~^idt<d ■ t^rh)) = M(7)"^"r(7) 

(7) 

= ^(7) ^ (u • fr(7)) = u - - N(7)Vs7 on (0, a) x E, 

which is the hidden evolution equation for 7. 

To keep this section short, we refrain from giving all details of the involved 
differential geometric computations here. For a coincise introduction of this topic 
we refer to [5], see also [5J Section 2]. 
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2.3. The Linearized Operators. Our analysis of problem ([5]) relies on maxi- 
mal Lp-regularity results for suitable linearizations. To obtain such a linearization 
w. r. t. reference functions 



u* e X„(a) 



c* € X,(a) := 



n Lp{{0, a), H^in \ E)") ■ [$]so = J ' 

e H^iiO, a), Lp(ext(E))) £ X,(a), 

n Lp{{0, a), 7l2(ext(S))) ' d,(j) = 



c* e X,(a) ffi((0, a), n Lp((0, a), 

first observe that 

M"ij) = dt + {u- dtQ) ■ 5(7) = dt + {u* ■ V) - M(?2, 7 | m*) • V, 
Gij) = (1 - G(7))V, 2?(7) = (1 - G(7))V • - div - 0(7) : V, 
C{j) = I?(7)5(7) = A - A(7) • V - L(7) : V^, 
where the multiplication operators are given as 

M(u, 7 I u*) = 9(9 - (u - u*) + D(7)(u - dtO), 

G(7) = (ve)-Hve-i), d(7) = g(7)t, 

A(7) = -(AO) o e, L(7) - D(7) - D(7)G(7) + G(7). 

Unfortunately, the transformation of the operators that act on functions defined on 
gr(r( • )) is somehow more involved. First note that a straight forward computation 
shows that 

5rh) = Pr(7)N(7)VE =: (1 - Gs(7))Vs. 

Indeed, we have 

Pr(7) = 1 - t'r(7) «) i^ril) = 1 - ^^{lfi''T. - N(7)Vs7) ® i'^T. - N(7)Vs7) 
= 1 - (1 - n(7)) (j/s - N(7)Vs7) «) {fs - N(7)Vs7) =: ft - P(7), 
where we have set 

|N(7)V.7P 



l + |N(7)Vs7p 
and, therefore, 

Gs(7) = (l-N(7)) + (l-Pr(7))N(7). 
To transform the material derivative for functions defined on gr(r( • )) we assume 

(j)r{t, ■) --Tit) — >R, 0<t<a 

to be a generic family of scalar valued functions. Setting 

hit, ■ ) := 0r(t, • ) o e(i, • ) : S — > R, 0<t<a 
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we obtain for fixed < t < a and a; G S the relation 

= pr0r)(t, e(t, x)) 
d 



ds 



s=t 



= ^'?^r(s, 0(s, a;)) 



+ ^^{<f>ris, X"(s 1 1, eit, x))) - 0r(s, e(s, a;))) 



= dt4>T.{t, x) + (us • Qt{i)) 4'T.{t, x), 

where u^, = P^u denotes the tangential velocity on the reference interface E and 
X"( • I i, x) denotes the characteristic curve of a particle, which is advected by the 
velocity field m, as defined in the introduction, i. e. 

ns\t,x) = u{s,x''{s\t,x)), \s-t\<e, x''{t\t,x)^x. 

Indeed, we have 

X"(s 1 x)) - ©(^, x) = (s - t) "(s, X"(s 1 e(t, x))) + o{s - t) 

=: (s — t) hi{s 1 1, x) 

e(s, x) - Q{t, x) = {s - t) 7(g. x)--i{t, x) ^^^^^^ 



s - t 
(s - t) ho{s 1 1, x) 



with 



lini^/ii(s \ t, x) — ho(s \ t, x)j = u{t, x) — dtj{t, x)i'^{x) 

= PT.{x)u{t, x) = U^{t, x) 

due to the kinematic condition ([7]) and, hence, 

^ (</)r(s, 1 1, e{t, x))) - 0r(s, e(s, x))^ 

^ Iim ('^r(s,x''(^|t,e(t,x)))-0r(t,e(t,a:)))-(0r(^.e(s,x))-0r(t,e(t.a^))) 



^ lim + ^1 3:))-0r(a. e(t, x) + (a-t) hojs, t, x)) 



= lim 



{hi{s \t,x)~ ho{s 1 1, x)) • (Vr0r)(<, e(t, x)) 



= (us ■t^r(7))0s(i, a;) 
Hence, the transformed material derivative reads 

with 

Ms(u, 7 I u*) = Ds(7)us - ("s - u'^), De(7) = Gs(7). 
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Finally, we need to linearize the curvature 



tr(N(7)(VsN(7)Vs7) + is) 

- Ai(7)^N(7)2Vs7- (VeN(7)Ve7)N(7)Vs7, 



which becomes 

k'(0) = trL| + As, 

cf. [S] and Section 2]. 

2.4. The Principal Linear Part. Based on the above preparations the principal 
linearization of problem ([5]) reads as follows. Note that here we drop the bars that 
have been introduced in the previous subsections for convenience. 

pdfU + (u* ■ V)m - TjAu + Vp = pFu{u, p, 7 | u*) in (0, a) x \ E, 

divu = G{u) in (0, a) x f7 \ E, 

dtc + {u* ■ V)c - dAc = 7, c | u*) in (0, a) x ext(E), 

Hs = on (0, a) X S, 
-PsI?7(Vu + VuT)]s//s 

- a'(c£)Vscs = He, 7, cs |c£) on (0, a) x S, 

(8) 

+ IpIs - fT(c^)Ascs - 7, CE I 4) on (0, a) x S, 

9(7 - u • t^E + (u£ • Vs)7 = H^{u, 7 I M*) on (0, a) x S, 

a'([clE)[c]E-cs - H^{c\c*) on (0, a) x S, 

9tCE + (u£ • Vs)cs - drAECs = H^{c, ce | u*) on (0, a) x S, 

Nao = 0, d^c ==0 on (0, a) x c)ri, 

m(0) — uq in il, 7(0) = 7o on E, 

c(0) — Co in ext(S), 

where the non-linearities are given as 

pFu{u,p,j\u*) = {pM(u, 7|u*) • V-77A(7) • V-77L(7) : V2}u 
+ {G(7)V}p, 
G(«) = {D(7): Vjii, 
7, c I m*) = { M(u, 7 I u*) • V + rfA(7) • V + dL(7) : } c 
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as well as 

Hliu, Iph, 7, c^\u*) = PsI?7(G(7)V^/ + G(7)VwT)]5.N(7)Vs7 

^ Ps[^(Vm + VuT)15.N(7)Vs7 
+ Ms N(7)Vs7 - ^sI?7(G(7)V?i + G{j)\/u^)]s 
+ M7)-Ml-M7))^'(4)Vscs 
+ M7)-M^'(cs)-rT'(4))Vscs 

- ^(7)"^ '^'(cs)Ge(7)Vscs 

- o-(cs)Kr(7)N(7)Vs7i 

7, cs I u*) = Iry(G(7)Vu + Gij)Vu^)j^ N(7)Vs7 • 

- |77(Vu + VuT)]s N(7)Vs7 • J^s 

- I?7(G(7)Vu + G(7)VuT)]5, t,^ . „^ 

+ KlY^ cr(cE)('^r(7) - As7)^'r(7) • vt. 

- ^J.{jy^ cr'(cs)GE(7)VsCs • I^E 

+ Kiy^ (f^(cE) - (T(4))Kr(7)i'r(7) • z^s, 
7 h*) = { (1 - N(7))us • Vs - (^/E - • Vs } 7, 
i7,"(c|c*) = a'([c*]E)[c]s-a([c]E), 
H^{c, CE I u*) = dd^c - cs X'r(7)u 

+ { Me(w, 7 I u*) • Ve - dAs(7) • Ve - rfLE(7) } cs- 

As mentioned above, our analysis of problem ([S]) - or equivalently problem ^ - 
relies on maximal Lp-regularity results for suitable linearizations, which here are 
given by the linear left-hand side of problem ([S]). These shall be established now. 
To be precise, we will prove the existence of a unique maximal regular solution 

(u, p, 7, c, ce) G X(a) := Xt,(a) x Xp(a) x X^{a) x Xc(a) x Xs(a) 

with 

Xp(a) := {qe Lp((0, a), : Ms G H;;(a) } , 

X^(a) := VFp'-^/'^((0, a), Lp(E)) n ^^((O, a), W^-'^^m 

n Lp((o, a), w^-'/'{j:)) 
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and X„(a), Xc(a) and Xs(a) defined as above to the linear problem 
pdtu + (u* • V)ii — ryAii + Vp = pfu in (0, a 

div u = g in (0, a 

dtc + {u* ■ V)c - dAc ^ fc in (0, a 

Hs = on (0, a 

-PsH^u + ViiT)]s j/s - (t'(4)Vscs = /ij; on (0, a 

(9) -l77(Vu + V7/T)]5.;/5..;/5. + |p]j,_a(c*)As7 = on(0, a 

9t7 — u ■ i^s + (u^ ■ Vs)7 = hj on (0, a 

a'([c*]E)[c]s -cs - /i? on (0, a 

9tcs + (us • Vs)cE - rfrAscs ft-^ on (0, a 

[u]on = 0, d^c = on (0, a 



X S1\S, 

X r2\ s, 

X ext(E), 
X S, 
X E, 
X E, 
X S, 
X S, 
X S, 

X 



u(0) = uq in 17, 7(0) = 7o on E, c(0) = cq in ext(S), 

provided the data satisfy the necessary regularity and compatibility conditions. 
Here the necessary regularity conditions are given as 



fu 


e F„(a) 


:= Lp((0, a), 


9 


e G(a) 


:= i7pH(0, a), oH-^^)) n ^^^'((0, a), Lpift)) 
n Lp((0, a), 


fc e F,(a) 


Lp((0, a), Lp(ext(I]))), 


hi 


e H-(a) 


iypi/2-i/2p^^Q^ ^^^^^ TS))nLp((0, a), W^-'/"" 


h'l 


e H-(a) 


M^p^/'-^/'P((0, a), ip(S))nLp((0, a), W^^'/^(S)), 




e H..y(a) 


iyp^^^/'^((0, a), nLp((0, a), W-p'-'/P(S)), 


h^ 




Xs(a), 


hi 




Lp((0, a), Lp(E)), 


Uo 


e T„ 


:= { $ e W^-^^^in \ E)" : = 0, ^an = } , 


70 






Co 




M^p'-'/P(ext(E)) : e M^p'-'/''(E), 9,0 = 



Note that all these spaces arise based on trace theorems for the solution space X(a). 
The first regularity assertion for the right-hand side g of the divergence constraint 
is a result of partial integration. Here oH~^{^l) := Hp,{iiy with l/p+ l/p' = 1. 
With the abbreviations 

Y(a) := F„(a) x G(a) x Fc(a) x H;(a) x H^(a) x H.^(a) x H^(a) x H^(a), 

T T„ X T..^ X Tc 
our main result concerning ([9]) now reads as follows. 
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Theorem 2.1. Let a > 0, let ^ M" be a bounded domain with boundary of class 
C^~, let S e M'H'^i^) be of class C^^ and let p > n + 2. Let p±, ?y±, d, dr > 
and let a, a £ C'^^(M+, IR+) with a' > 0. Moreover, let 

u*eXu{a), c* eXcia), c*^eXs{a). 

Then the linear problem admits a unique maximal regular solution 

{u, p, 7, c, cs) e X(a), 

if and only if the data satisfy the regularity conditions 

ifu, 9, fc, K, hi, h^, /i^, hi) e Y(a), K, 70, Co) G T 

and the compatibility conditions 

divuo = g(0), -FeI?7(Vuo + V4)ls - a'(4(0))VsSo = /^^(O) 

So = a'([c*(0)]s)Ns"/i?(0). 
Moreover, the solutions depend continuously on the data. 

Proof. We start by determining cs G Xs(a) as the unique maximal regular solution 
to the parabolic problem 

dtCY. + • Vs)cE - rfpAscs = /i^ on (0, a) x S, 
cs(0) = So on S, 

which then enables us to obtain c G Xc(a) as the unique maximal regular solution 
to the parabolic problem 

dtc + {u* ■ V)c - dAc = /c in (0, a) x ext(i;), 

a'([c*]s)[c]s = cs + /i? on (0, a) x E, 

9,yC = on (0, a) x dVl, 

c(0) Co in ext(i;). 

Note that at this point we require a' > in order to force the boundary condition 
to satisfy the Lopatinskii-Shapiro condition. Now, the remaining system for the 
unknowns 

{u, p, 7) G X„(a) X Xp(a) x X-y(a) 

forms a two-phase Navier-Stokes problem with time and space dependent surface 
tension cr(c|^) and we may invoke [BJ Theorem 3.1 & Corollary 3.4] to complete the 
proof. □ 
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3. Local Well-Posedness 

3.1. The Non-Linear Problem. With the maximal Lp-regularity result, Theo- 
rem [2lll at hand we require only one more ingredient to prove the well-posedness 
assertion of Theorem ll.il Of course, we also need suitable estimates for the non- 
linear right-hand sides in ([8|). To shorten the notation we denote by 

L(-|u*, c*, :X(a) — >Y{a), resp. N{- \u* , c* , c*^) : X{a) — >Y{a) 

the linear operator induced by the left-hand side resp. the non-linear operator in- 
duced by the right-hand side of ([8]) without the inital condition. Thus, the problem 
([5]) is equivalent to 

L{u, p, 7, c, cs I u*, c*, c^) = N{u, p, 7, c, cs | u* , c*, c^), 

{u{0), 7(0), c(0)) = (uo, 70, Co) 

and the operator L{ ■ \ u*, c*, c|.) is bounded, since the derived regularity conditions 
on the data are necessary. Moreover, the operator 

oL( • I u*, c*, Cs) : oX(a) — > oY{a), 

which denotes the restriction of the operator L{- \u* , c* , c'^) to the closed sub- 
space oK{a) of X(a) with vanishing initial values, is an isomorphism thanks to 
Theorem 12.11 Here we denote by Y(a) the closed subspace of Y(a) with vanish- 
ing initial values and note that the corresponding compatibility conditions on the 
data are necessarily satisfied in the setting of vanishing initial values. Now, the 
non-linear operator N{- \u*, c* , c^) has the following mapping properties as has 
been shown in [71 Proposition 4.1] for the two-phase Navier-Stokes equations with 
surface tension in the case where S = R"~^ and in [SJ Proposition 4.2] in the case 
of a general geometry as considered here. 

Proposition 3.1. Let a > 0, let Q R" he a hounded domain with houndary 0/ 
class , let E g M'H'^{Vt) he of class and Ictp > n+2. Let p±, r]±, d, dr > 
and let a, a e C^^(M+, M+) with a' > 0. Moreover, let 

u* G Xi,(a), c* G Xc(a), G Xs(a). 

Then 

N{-\u\c\cl)eC\X{a),Y{a)) 
and the Frechet derivative DN{u, p, 7, c, cs | m*, c*, d^) satisfies the estimate 

M-^\\DN{u, p, 7, C, Cs I U*, C*, c£)llB(oX(a),oY(a)) 
< II U - U*|lBC((0,a),BCi) + 11"^ ^ ^*\\BC{(a,a),BC^) 
+ I|CE - Cslli3C((0,a),BCi) + IK"' 7, C, Cs)||x(a) 
+ (Il7llx^(a) + l|Vs7llBC((0,a),i3Ci))ll("' 7, C, CE)||x(a) 
+ ^(Il7llx-,(a)' llVs7lli3C((0,a).BCi))ll7llx-,(a) 
+ (3(11 Vs7||bc'((0, S), BCi)) ll'^S7lli3C((0, a), SCI) 
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for all a G (0, a] and all (m, p, 7, c, cs) G X(a) wzt/i some constant M = M{a) > 0. 
Here P and Q denote fixed polynomials. Moreover, if a, cr G C"^(R4_, then 

N{-\u*,c*,c*^)eC^{X{a),Y{a)). 

The proofs given in the above mentioned sources carry over to the case considered 
here. The main arguments are the polynomial structure of N{ ■ \ u* , c* , c'^), which 
is equal to that considered in [6^ , and the embeddings 

X„(a), X^(a), Xc{a) BUC{{0, a), BUC^), 

which are available thanks to p > n+2 and which imply the involved function spaces 
to be Banach algebras. With this preparations at hand the proof of Theorem ll.il 
may be carried out as follows. 

3.2. Proof of Theorem II.IL Step 1. Suppose the assumptions of Theorem ll.il 

are satisfied and 

uo G W^-'/Pin \ To), To G CO G w^-'/pin+{o),R+) 

that satisfy the stated regularity and compatibility conditions are given. According 
to the considerations in Section [2] we may approximate To for every given 6 > 
with a real analytic hypersurface E G M7i^{il) such that 

distA4W2(o)(I], To) < S 

and Fq is parametrized over E by a function 70 G Wp~^^P (T,). Therefore, it is 
sufficient to consider the transformed problem ([8]), where we drop the bars again 
for convenience. 

Now, we denote by u G X,j(a) an arbitrary extension of uq, which may e.g. be 
obtained as a solution to a suitable parabolic problem. Moreover, we denote by 
G EI"(a) an arbitrary extension of a'(0)[co]E, which may also be obtained as a 
solution to a suitable parabolic problem on S. Finally, we denote by g G G(a) resp. 
h'^ G El^(a) extentions of 

divMo resp. - -Ps [7?(Vuo + Vj)]s i^s, 

which exist thanks to [SI Proposition 4.1], and solve 

{u\ p\ 7*, c\ 4) G X(a), (u*(0), 7*(0), c*(0)) = K, 7o, cq) 

L{u\ p*, 7*, c*, 4 I u, 0, 0) = (/, g, 0, hi, 0, 0, K, 0) 

to obtain a reference solution [u* , c* , d^) G X„(a) x Xc(a) x Xs(a). 

3.3. Proof of Theorem ll.il Step 2. To prepare a fixed point argument we split 
the desired solution z ~ (u, p, 7, c, cs) G X(a) as 

z = z* + z, z G oX(a), 

where z* = [u* , p* , 7*, c*, c^) G X(a) denotes the reference solution obtained in 
the first step and a G (0, a] will be chosen below. Problem ([5]) may then be rewritten 
as 

oL(z \u*,c*, c£) = N{z + z* I u*, c*, 4) - L{z* \ u* , c* , 4) =: K{z) 
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and, hence, the solution is given as the fixed point z = oL{- \ u*, c* , c^) ^K{z). 
Now, 

\\oL{-\u*, C*, c£)"^|ls(oY(a),oX(a)) < C*, aG (0, a] 

and due to K{0) =^ N{z* \ u*, c*, c'^) ~ L(z* \ u*, c*, c^) and Proposition 13.11 we 
may choose a G (0, a] and r > sufficiently small to ensure 

r _ 1 _ _ _ 

Thus, the contraction mapping principle applies and yields the unique solution 

z e o'K{a). 

4. Analyticity 

4.1. Proof of Theorem Step 3. In order to complete the proof of Theo- 
rem [TTT] it remains to show the analyticity of the solutions. However, this may be 
obtained by the well-known parameter-trick since the non-linear right hand side 
N{ ■ I u*, c*, c^) is real analytic, provided a, a e C"(R+, R+), cf Proposition 13. II 
We refrain from giving all the details here, since the proof is similar as that of [SJ 
Theorem 4.3] and [3 Theorem 6.3]; see also [3l|4]. 

5. Semiflow 

5.1. Proof of Theorem II. 2i That the solutions to ^ obtained by Theorem ll.il 

generate a local semiflow in the phase manifold Sp{fl) can be seen as follows. Given 

u„ e w^''/p{n \ Fo), Fo e w^-^/p, CO G w^-'/pin+io),R+) 

that satisfy the regularity and compatibility conditions as stated in Theorem ll.il 
we obtain a local solution on some time interval [0, a) with a > 0. However, 
these solutions belong to the maximal regularity class defined in Section [5] and, 
thus, admit a trace at t = a, which allows to reinvoke Theorem 11.11 to continue 
the solution on a larger time interval. Note that this also allows to choose a new 
reference manifold for the construction of the solution. Repeating this procedure 
we obtain a maximal time interval [0, a*), which may be bounded by the fact that 

lim {u{t), F(t), c{t)) 

does not exist in Sp{fl) or by the fact that the interface undergoes a change in its 
topology or touches the outer boundary 9r2, in which case the model is no longer 
suitable to describe the situation. 
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